Thermal width and gluo-dissociation of quarkonium in 

pNRQCD 



Nora Brambilla a , Miguel Angel Escobedo a , Jacopo Ghiglieri a ' b , Antonio Vairo a 

" Physik-Department, Technische Universitat Miinchen, James-Franck-Str. 1, 85748 

Garching, Germany 

b Excellence Cluster Universe, Technische Universitat Miinchen, Boltzmannstr. 2, 85748, 

Garching, Germany 



Abstract 

The thermal width of heavy-quarkonium bound states in a quark-gluon plasma 
has been recently derived in an effective field theory approach. Two phenomena 
contribute to the width: the Landau damping phenomenon and the break-up 
of a colour-singlet bound state into a colour-octet heavy quark-antiquark pair 
by absorption of a thermal gluon. In the paper, we investigate the relation be- 
tween the singlet-to-octet thermal break-up and the so-called gluo-dissociation, 
a mechanism for quarkonium dissociation widely used in phcnomenological ap- 
proaches. The gluo-dissociation thermal width is obtained by convoluting the 
gluon thermal distribution with the cross section of a gluon and a 15 quarko- 
nium state to a colour octet quark-antiquark state in vacuum, a cross section 
that at leading order, but neglecting colour-octet effects, was computed long 
ago by Bhanot and Peskin. We will, first, show that the effective field theory 
framework provides a natural derivation of the gluo-dissociation factorization 
formula at leading order, which is, indeed, the singlet-to-octet thermal break-up 
expression. Second, the singlet-to-octet thermal break-up expression will allow 
us to improve the Bhanot-Peskin cross section by including the contribution of 
the octet potential, which amounts to include final-state interactions between 
the heavy quark and antiquark. Finally, we will quantify the effects due to final- 
state interactions on the gluo-dissociation cross section and on the quarkonium 
thermal width. 

Keywords: Quarkonium, finite temperature, thermal width, gluo-dissociation, 
singlet-to-octet break-up 



1. Introduction 

Quarkonium suppression has been suggested long ago as a hard probe of 
the medium produced in heavy-ion collisions [1 . This hypothesis has been 
widely investigated, both theoretically and experimentally, in the past 25 years 
[21 |3]. The early theoretical arguments were based on the expectation that 
above the deconfinement temperature the linear, confining part of the heavy 
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quark-antiquark (QQ) potential would vanish and the Coulomb part at short 
distances would be replaced by a screened Yukawa (or Debye) potential that can 
support only a limited number of bound states. Since the screening (Debye) 
mass depends on the temperature, heavy-quarkonium states were thought to 
provide a thermometer of the medium. 

In the last few years, significant progress has been made in deriving the 
heavy QQ potential from QCD in a rigorous and systematic way. The real-time 
static potential was first calculated for large temperatures, T ^> 1/r > mjj, 
where r is the quark-antiquark distance and mo is the Debye mass, in [3HZ] - 
For a wider range of temperatures, an effective field theory (EFT) study of non- 
relativistic bound states at finite temperature has been carried out for QED in 
[5HTU] and for QCD in [TTMl3| . Most importantly, in the same framework, also 
the quarkonium thermal width has been calculated. Two mechanisms, at least, 
have been identified as responsible for it: the Landau damping phenomenon |3] 
and the singlet-to- octet thermal break-up [11] . In the former, the virtual gluons 
that are exchanged between the QQ pair scatter off the light constituents of the 
medium, whereas in the latter the colour-singlet bound state absorbs a gluon 
from the medium and turns into a colour-octet state. 

There exists, however, a large literature where the quarkonium behaviour in 
a medium is studied on a phenomenological basis (see e.g. (T3Hl6| ). The thermal 
decay width is obtained by convoluting scattering cross sections computed at 
T = with thermal distributions for the incoming light partons. At least 
two scattering processes have been considered: gluo- dissociation and quasi-free 
dissociation. It is then natural to ask, if and to what extent gluo-dissociation and 
quasi-free dissociation agree with the singlet-to-octet break-up and the Landau 
damping widths derived from EFTs. 

In this letter, we will deal with gluo-dissociation [T7j[18] (see [T5H36] for some 
recent literature and [TSlIiniEZ] for reviews). The mechanism of quasi- free disso- 
ciation 38] and its relation with the EFT framework will be dealt with elsewhere 
[59"] . The process underlying gluo-dissociation is the same that gives rise to the 
singlet-to-octet thermal break-up width in the EFT framework. The cross sec- 
tion for gluon absorption by a colour-singlet 15 state was computed in 1979 by 
Bhanot and Peskin (BP) in Refs. [301111]: where the contribution from the final- 
state interactions was neglected by considering the large number of color, N c , 
limit. In the present work, we will (i) prove at leading order the factorization 
formula for quarkonium gluo-dissociation and establish under which conditions 
it holds, hence, show that the gluo-dissociation thermal width coincides with the 
singlet-to-octet thermal break-up width, (ii) improve the gluo-dissociation cross 
section by including the contribution of the octet potential, which amounts to 
include QQ final-state interactions, and, finally, (Hi) assess quantitatively the 
effects due to final-state interactions on the gluo-dissociation cross section and 
on the thermal width. 

The letter is organized in the following way. In the next section, we will recall 
some basics on gluo-dissociation. In Sec. [3j we will derive the gluo-dissociation 
factorization formula in an EFT framework, indeed showing that it coincides 



2 



with the singlet-to-octet thermal break-up expression, first derived in [TJ]. We 
will also show that the gluo-dissociation cross section agrees with the BP cross 
section in the large N c limit. In Sec. [4] we will improve the BP cross section by 
including final-state interactions between the QQ pair in a colour-octet state. 
Finally, in Sec. [5j we will draw some conclusions. The results presented here 
are also part of the Ph.D. thesis [32] • 

2. Gluo-dissociation 

Gluo-dissociation stands for the process g+$(lS) — > {QQ)s, where a quarko- 
nium IS state, $(15), absorbs a gluon and becomes an unbound QQ pair in 
a colour-octet state, (QQ)s- hi the literature, it has been assumed that con- 
volving the in vacuum gluo-dissociation cross section, a±s, with the thermal 
distribution of the gluons provides the gluo-dissociation thermal width, risj^] 
Specifically, one writes (see for instance Eq. (23) of [T5] ) 

Tis=/ 7^n B (|q|)a lff (|q|), (1) 

J\1\>\E 1S \ ( 271 T 

where ns(x) = {e x / T — l)^ 1 is the Bose-Einstein distribution, E\s the binding 
energy of the quarkonium 15 state, and we have assumed the bound state and 
the bath to be at rest. 

In [301 |3I] , the gluo-dissociation cross section was calculated at leading or- 
der, under the following assumptions: (1) the quarkonium 15 state, $(15), 
is Coulombic; (2) in an operator product expansion framework, the gluon- 
quarkonium interaction is taken at leading order, which corresponds to a chro- 
moelectric dipole interaction; (3) the (repulsive) octet potential is neglected, 
which is tantamount to neglecting final-state interactions. Theoretically the 
last assumption may be realized by taking the large-7V c limit. In this limit, 
the colour-singlet Coulomb potential, i.e. the potential between a QQ pair 
in a colour-singlet configuration, which is (at leading order) v} ^ = —CfO-s/t 
with C F = {Nj - 1)/(2JV C ), becomes v}% = -N c a s /(2r), whereas the colour- 
octet Coulomb potential, i.e. the potential between a QQ pair in a colour- 
octet configuration, which is Vo — a s /(2N c r), vanishes. The large-iV c limit 
also modifies the Bohr radius from o,q = 2/(mCfa s ) to obp = 4/(miV c a s ) = 
4/(3ma s ) and the absolute value of the binding energy (at leading order) from 
| i?i | = mCpOt^/4 to £i,bp = l/( ma lp) = 9m«g/16. The Bhanot-Peskin gluo- 
dissociation cross section, cris,Bp(|q|), as a function of the gluon momentum q, 
then reads 

^,BP(|q|) - — g m • (2) 



1 There exist papers, such as [231 OS], where finite-temperature effects in the cross section 
are included by considering colour-singlet wavefunctions derived from potential models. On 
the thermal distribution side, hydrodynamical and anisotropic effects have been considered in 
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We note that, in the above formula, an overall colour factor Cp = 4/3 has 
been kept unexpanded when performing the large N c limit. We call Tis^p the 
corresponding gluo-dissociation thermal width. 

The cross section |2]) is averaged over the 2 polarizations and the 8 colours 
of the initial gluon. Hence, when inserted in the decay width formula, Eq. 
the Bhanot-Peskin cross section should be multiplied by a factor 16: <ris(|q|) « 
16 x ois,Bp(|q|)- The factor is explicitly included in Eq. (4) of [25] and also the 
authors of [33] multiply the BP cross section by 16 when using Eq. fllr] 



3. Effective field theory 

The EFT approach is based on the hierarchies of non-relativistic and thermal 
scales typical of quarkonium in a quark-gluon plasma. The hierarchy of non- 
relativistic scales follows from the fact that the heavy quark has a velocity in the 
centre-of-mass frame that is v 1 ; the hierarchy is then m 3> mv 3> mv 2 , where 
m is the heavy-quark mass, mv is the scale of the typical momentum transfer in 
or inverse radius of the bound state and mv 2 is the scale of the typical energy. 
For a Coulombic bound state, such as the bottomonium ground state likely is, it 
holds that mv ~ ma s 3> Aqcd and also that E ~ ma 2 > Aqcd- The hierarchy 
of thermal scales is T 3> mo, where T is the temperature of the quark-gluon 
plasma and mo the Debye mass|^] 

The relative size of non-relativistic and thermal scales depends on the medium 
and on the quarkonium state. In the following, we will adopt the hierarchy 

m ^> mv ~ ma s ^> T ~ mv ~ ma 2 ^> mp, Aqcd , (3) 

which is meant to include also the regions mv ~ ma s 3> T > mv ~ ma 2 and 
mv 2 ~ ma 2 >T^> mo, Aqcd, although this last one is of less phenomeno- 
logical impact for the thermal width is exponentially suppressed |llj . The hi- 
erarchy ([3]) is more general than the one analyzed in |12j . where we required 
mv ~ ma s T mv 2 ~ ma 2 . It was argued in [43 that this hierarchy may 
be the relevant one for T(15) produced in heavy- ion collisions at the LHCnln 
|12j . the spectrum and width of quarkonia were computed up to order maf^ 

EFTs suitable to describe quarkonium in a medium are constructed by sub- 
sequently integrating out high-energy scales in ([3]). Integrating out modes that 
scale like m and ma s leads respectively to non-relativistic QCD (NRQCD) 
[46j 147] and potential non-relativistic QCD (pNRQCD) [48] [49]. Since the 



2 Private communications from Xingbo Zhao are acknowledged. 

3 The right temperature scale is rather irT, or multiples thereof, than T. This is, however, 
relevant, only when quantifying the different energy scales in the system; hence, we will drop 
the factor n from qualitative considerations. Effects due to the magnetic mass are suppressed 
and do not contribute to the considered accuracy. 

4 We refer to 1441 for the most recent CMS measurements on the suppression of the T 
family. A phcnomenological analysis of the data that includes the effects of the Landau- 
damping width can be found in [45 j ■ 



4 



temperature is much smaller than both m and ma s , it can be set to zero in 
the matching and both Lagrangians are the same as at zero temperature. The 
pNRQCD Lagrangian, in particular, reads 



£ P NRQCD = -\F; v F a ^ v + J2q z lpq t 

i=l 

d 3 r Tr {S f [ido -h.]8 + O f [iD - h a ] O 

+V A (O t r-.gES + H.c.) + {r • 5 E , 0} + . . . }. (4) 

The fields S = S lc/VN~c and = O a T a j \jTp are the QQ colour-singlet and 
colour-octet fields respectively, n/ is the number of light quarks, qi, Tp — 1/2, 
E is the chromoelectric field, iD O = idoO — gA O + OgA and H.c. stands 
for Hermitian conjugate. The trace is over colour and spin indices. Gluon fields 
depend only on the centre-of-mass coordinate and on time; this is a consequence 
of having multipole expanded the gluon fields in the quark-antiquark relative 
distance r. The dots in the last line stand for higher-order terms in r and 1/m. 
We note that, as in the BP approach, the leading gluon-quarkonium interaction 
is a chromoelectric dipole interaction. 

The dependence on the scales m and ma s is encoded in the Wilson coeffi- 
cients; Va and Vb are at leading order Va = Vb = 1, whereas the singlet and 
octet Hamiltonians have the form (p = — «V r ) 

n 2 V {1) V {2) 

m m m z 

The dots stand for higher-order terms in the 1/m expansion. The first two terms 
in the right-hand side, which are the kinetic energy and the static potential, 
constitute the leading-order Hamiltonian. The singlet and octet leading-order 
Hamiltonians explicitly read 

hV> = £-c P * *g», = p! + i*. (6) 

m r m 2Jy c r 

The spectrum of hs is made by the (QCD) Bohr levels E n = —mCpal /(An 2 ), 
whereas the octet potential is repulsive and does not support bound states 
but a continuum of scattering states. Note that, in the non-relativistic EFT 
power counting, both the kinetic energy and the static potential scale like ma 2 . 
Therefore, neglecting the octet potential, as done in the BP calculation, is a 
sensible approximation only in the large N c limit. 

We now set out to compute the singlet-to-octet break-up amplitude and 
width within pNRQCD. We refer to [TTJ El for details regarding pNRQCD at 
finite temperature in the real-time formalism]^] The singlet-to-octet break-up 



5 Since we do not require here T 3> mOj 2 , a difference with 1121 is that we avoid integrating 
out T and constructing the intermediate EFT called pNRQCDuTL [50] . 
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Figure 1: The leading heavy-quarkonium self-energy diagram in pNRQCD. The single line is 
a singlet QQ propagator, the double line an octet QQ propagator, the curly line a gluon and 
the vertices are chromoelectric dipole vertices. The imaginary part, obtained by cutting this 
diagram, gives the singlet-to-octet break-up thermal width. 



thermal width is given by the imaginary part of Fig. [T] The amplitude £(£7) of 
that diagram reads 



£(£) = -ig 2 C F -r* 



d 4 k 



(2tt) 4 E - -k + i V 
x [kl2nn B (\k \)6{kl-k 2 )] r\ 



(7) 



where E is the energy of the incoming QQ singlet and the expression in square 
brackets is the thermal part of the chromoelectric correlator. Since T ~ E 3> 
mr>i bare propagators can be used. The number of dimensions has been set to 
4, the integral being convergent. Evaluating it yields for the imaginary part 



ImS(B) = 



9 2 C f 
6ir 



The singlet-to-octet break-up width for the 15 state then reads 



r ls = -2(LS|Im | IS) = 



9 2 C F 
3tt 



E-hM "n B {\E-h^\y. (8) 
ite then reads 

(ISy E 1 -h^\ 3 n B (\E 1 -h^\)r i \lS), 

(9) 

where (rjlS") = l/(v / 7ra 3 ^ 2 ) exp(— r/ao) is the Coulomb 15 wavefunction. The 
corresponding expression in the static limit was obtained in |llj and in QED, 
for the hydrogen atom, in |S]. 

The difficulty in the evaluation of Eq. ^ lies in the Bose-Einstein dis- 
tribution and in its nontrivial dependence on ■ In Ref. |12j . having in- 
stead assumed T 3> ma 2 , we could expand the Bose-Einstein distribution as 



n B (|£i - h, 
ma b s E 1 /T 



(0)| 



rif 



T/\E X - 
Tal 



1/2 + . . ., obtaining up to corrections of order 



(4CJ. + iC 2 F N c + C F N 2 



24 



WC 3 F 



20C 2 F N C 



+ 8C F N 2 + N?) 



(10) 



The terms in the first line are the leading ones and are linear in the temperature]^] 



The linear behaviour of the thermal width has been recently investigated and found 
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Without expanding the Bose-Einstein distribution, the matrix element in 
Eq. (9) can be evaluated analogously to how the QCD Bethe logarithms, which 
contribute to the quarkonium spectrum at order ma^, have been dealt with 
in [551 , i- e - by inserting a complete set of octet states. Octet states are 
labeled by their energy and angular momentum quantum numbers and obey 
h °\pU z ) = (p 2 /m)\pll z ). It is convenient to introduce an arbitrary unit vector 
p and define a state |pZ) = (An/p) \pllz){Hz\p)-i where (ll z \p) = (p)* 



spherical harmonics. A suitable normalization of the states |pZ) is 



E 



d 3 p 

W) 



3N x|pZ)<pZ|y) =6 3 (x-y). 



is a 



(11) 



Inserting (11) into ^ gives 

(]q|) 



is 



|q|>|Ei| 



3tt 



|p|=V m (kl+Ei) 



,(12) 



where we have used that 



E 1 



n B (\E v 



is analytic in ho and made 



explicit that (lS\r projects on a I = 1 state. Equation (12) provides, at leading 
order in the EFT power counting, the decay width associated to the quarkonium 
singlet-to-octet thermal break-up, which is the dominant decay process in the 
situation (f3l). Quarkonium singlet-to-octet break up describes, at the order we 
are working, the same process of quarkonium gluo-dissociation, and the thermal 
decay width ( [l2| ) may be identified with the gluo-dissociation thermal width. 
In the following, we will show that Eq. ( 12 ) satisfies, indeed, the properties of 
the gluo-dissociation width presented in Sec. 2l 

First, since the thermal decay width (fl2) is expressed as a convolution of 
the gluon Bose-Einstein distribution and a function of the gluon momentum, 
this proves, at leading order, the factorization formula ([I]). It also allows the 
identification 



°is(|q|) 



3tt 



— l(lS|r| P l)| 2 



|p|=Vm(|q|+Bi) ' 



(13) 



for |q| > \E\\. We note that, in an alternative derivation, we could have made 
use of cutting rules at finite temperature [54j for the imaginary part of the 
diagram in Fig. [T] This would have led again to the factorization formula ([I]) 
and to the identification of <r\s with a T = cross section. However, the 
factorization formula ([T]) is not expected to hold at higher orders. 



consistent with lattice data in I51| . 
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Finally, in order to reproduce the BP gluo-dissociation cross section ^ from 
(13), we evaluate the dipole matrix element squared | (lS'|r|pl) | 2 in the absence 
of the octet potential. This is tantamount to using plane waves for the octet 
wave functions: ^^(x|pZ) (pZ|y) = e lp '^ x_y \ The matrix element squared then 
i 

becomes the square of the derivative of the momentum-space wavefunction and 
reads 

|<L%| P 1)| 2 |V p (p|LS)| 2 = ( f^ff )6 . (14) 



Plugging this into Eq. (13) and replacing E\ — > — £i,bp, Cf — 4/3, we obtain 

h w c ->oo ic 2 9 7ra s 4% (|q| - £i,bp) 3/2 1C n ,., ^ 

o-is q — -> 16— ^ = 16cti5,bp q , 15 

9 m |q| a 

r iS f 7 ^3« B (|q|)16<7 ls ,Bp(|q|)=r 1 s,Bp. (16) 

Equation (15) reproduces the BP cross section as already discussed, the 
factor 16 accounts for the 2 polarizations and 8 colours of the gluon. Therefore, 
the EFT computation in the large- N c limit leads naturally to the BP factoriza- 
tion and cross section formulas for quarkonium gluo-dissociation in a medium. 
We now set out to include the octet potential in the calculation, thereby also 
quantifying the approximation introduced by neglecting it. 



4. Colour-octet effects 



The calculation of the dipole matrix element squared |(15|r|pl)| , when 
octet potential contributions are included, is more involved, and requires the 
explicit integration over the continuum octet wavef unctions. Coulombic wave- 
functions in the continuum region \pll z ) can be found in [55]; I = 1 octet wave- 
functions |pl) can be found in [52j[53]. After correcting some typos, they read 



r pi = e 



- p i(ir/2-«i), 



27rp • r 



P 1 



1?,\P\ 2 



~TTpT - i 



i| P |H 



^ a |p| (e»oi 

x 1 F 1 (2 + ip/(a |p|);4;-i2|p||r|), (17) 



where iF\ is the confluent hypergeometric function, Si is the I — 1 Coulomb 
phase and p = l/(iVj — 1). The matrix element squared is then 

2 5127r 2 p(p + 2) 2 ag|p| (l + ^) e 4k-^(a | P |) 
l<lS|r|pl)| 2 = , 2jp \ °' Pl = • (18) 



e =5lpT-l (l + a 2|p|2) L 
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It is easily seen that the N c — > oo (p — > 0) limit of this equation gives back 
Eq. (14 1. Plugging the matrix element into Eq. (13 1 yields 



°is(|q|) 



3 TO|q| 



;t(iqi) 2 +p 2 : 



exp (ntqir arctan WlqD) 



e*(iqi) — 1 



(!9 



where t(\q\) = v /jq|7jJ3ij - 1. The limit N c -> oo (p -> 0) gives back Eq. (ITsJ 




Figure 2: In the main plot, the full cross section <ri s , given in Eq. \19\ , is plotted in continuous 
black. The BP cross section, 16<tis bP! Eq. \15\ , with the substitution ei^ap — > \Ei\, is 
plotted in dashed red. The inset plot shows the ratio o"is/(16 crigBp)- The horizontal dot- 
dashed line is the asymptotic limit (17/16) 2 , which is reached from above. The dotted blue 
curve in the main plot is the BP cross section without the substitution ei bp — > 



In order to estimate the approximation introduced by ignoring the octet 
potential, in Fig. [2j we plot with a continuous black line the exact cross section 
( 19 ) and with a dashed red line the BP cross section ( 15 1 as functions of the gluon 



The 2S cross section reads 



°"2s(|q|) 



a B C F 



2 ld 7rV 



[2E 2 (2p 2 + 5p + 3) + [q[(p + 2)] 2 (t 2 (|q|) 2 + 4p 2 ) 



x ex P ( i /i in arctan (*2 (|q|)) 

,*2(|q|) 



where t 2 (|q|) == VM/I^I - 1. 



;t 2 (iqi) _ I 
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momentum |q|. For a meaningful comparison, for the latter we have performed 
the substitution e^ep — > \Ei\, which guarantees that the binding energies are 
the same in the two cases. The dotted blue line is the BP cross section (151 



without this substitution, i.e. in terms of e^BPi which shows a larger threshold 
and a smaller peak. In the inset plot, we show the ratio, <Jis/{^-Qo'is,bp), of the 
continuous black and dashed red curves: the horizontal dot-dashed line is the 
asymptotic value of (2 + p) 2 /4, which, for N c = 3, yields (17/16) 2 . 



£. 2 a s Cp I m 




Figure 3: The width T\g is shown as a continuous black line and the width Tis bp as a dashed 
red line. 

In Fig. |3J we plot the widths Y\s and Tis^p as a function of the tempera- 
ture. They have been obtained by numerical integration of Eqs. (12 1 and (16 1 



respectively. In the latter case, we perform the same substitution as above, i.e. 
£i.bp — ^ \E\\- Tig is the continuous black line, whereas Fi^bp is the dashed 
red line. One clearly sees how the full result T\s overpowers the old BP result 
in whole range T > \E\\, and how the two widths quickly reach an asymptotic 
linear regime for T |-Ei|, as predicted by the analytical result (10 1. For 
T < \Ei \ the two widths become exponentially small. 

The analytical, asymptotic expression of the full width Tis for T ^> |2£i| is 



in Eq. (10). In terms of p, it reads 



EfC F a s /r 



16 



(2 + P) 2 



T 



-(2 + p) 2 (3 + p) 



O 



m 

T 



(20) 



In Fig. |4j the width (20 1 for N c = 3 (p = 1/8), corresponding to the inclusion 
of the octet potential, is plotted as a dot-dashed black line. The corresponding 
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Figure 4: The width Fig is shown as a continuous black line and the corresponding analytical 
result for T> \Ei\, Eq- l |20| with p = 1/8, is plotted as a dot-dashed black line. Similarly, 
the width Tis bp is plotted as a dashed red line and the corresponding analytical result for 
T > |_Bi|, Eq. Q| with p = 0, as a dotted red line. 



plot of Tig, obtained from a numerical integration of Eq. ( 12 1, is the continuous 
black line. The width (20) for N c —> oo (p = 0), corresponding to the BP 
approximation of a vanishing octet potential, is plotted as a dotted red line. The 
corresponding plot of Tx^bPj obtained from a numerical integration of Eq. ( 16 ), 
is the dashed red line. Both Tig and F^ bp approach their asymptotic linear 
regimes starting from T s» 4|£i|. For T s» 4|_Ei|, Tis is still larger than its large 
T asymptotic value by about a factor 2.7. In Fig.[5j we plot the ratio Tis /^is,bp, 
which also shows the deviation from the asymptotic limit of (17/16) 2 « 1.13. 

Finally, we remark that our results are valid as long as ma s 3> T (see Eq. 
([3])), which guarantees that the interaction of the bound state with the thermal 
gluons can be treated as a chromoelectric dipole. Therefore the regions of the 
plots where T /\E\\ > l/a s ~ 5 are to be intended for illustration purposes only. 



5. Conclusions 

In this letter, we have shown that, under the scale hierarchy ([3|, the leading 
contribution to the thermal width of a quarkonium IS* state may be written as 
a convolution integral of the gluon distribution function and a cross section ais, 
see Eq. (JlJ . The underlying process is known as singlet-to-octet thermal break- 
up in the EFT literature and as quarkonium thermal gluo-dissociation in the 
phenomenological literature: a colour-singlet QQ state interacts with a gluon of 
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0.0 



T 

I s, 



Figure 5: Plot of the ratio Tis/^lS bp; the horizontal dot-dashed line is the asymptotic limit 
(17/16) 2 . 



the thermal bath and breaks up in an unbound colour-octet QQ pair. The cross 
section <7is can be identified with the in vacuum gluo-dissociation cross section. 



We have derived its explicit expression in Eq. (19). This expression includes, 
for the first time, the contribution of the colour-octet potential, i.e. the final- 
state interactions between the heavy quark and antiquark. We have shown that 
the gluo-dissociation cross section and decay width reduce to the well-known 



Bhanot-Peskin result if colour-octet effects are ignored, see Eqs. (151 and (16) 



Under the condition T 3> raa s 2 , the thermal decay width may be expanded to 



give back the singlet-to-octet break-up width calculated in [12], see Eq. (10). 

The EFT framework, in which the factorization formula and the cross-section 
expression have been derived, makes clear the region of validity and the accu- 
racy of the obtained results. They hold under the conditions ma s ^> T, so 
that the interaction of the bound state with the gluons of the medium can be 
described by a chromoelectric dipole interaction, mai? >mi) and T mo, so 
that, in first approximation, the thermal masses of the gluons can be neglected. 
The factorization holds at leading order in the EFT power counting. Beyond 
leading order, which includes contributions coming from hard thermal loop re- 
summed gluon propagators, the simple convolution formula will break down. 
The explicit expression of the gluo-dissociation cross section that we have pre- 
sented, as well as the old Bhanot-Peskin expression, follows from the assumption 
that the QQ pair is weakly coupled so that the colour-singlet bound state and 
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the colour-octet unbound quark-antiquark pair may be described in terms of 
Coulombic bound or scattering states respectively. This assumption is likely to 
hold only for the quarkonium ground state; for this reason, we have restricted 
our analysis to the dissociation of quarkonium IS" states, although an extension 
to quarkonium states with arbitrary quantum numbers would be straightforward 
(the cross section for an arbitrary state can be found making the substitutions 
Ei -> E n , (1S\ -> (nll z \ and |pl) -> J2i\P l ) in Ec l- @> see footnote [7] for 
the explicit expression of the 02s cross section and Fig. |6j for the T 2 s width). 
In view of this, we stress that phenomenological studies that make use of the 
Bhanot-Peskin gluo-dissociation formula, but fix the binding energy by some 
non-Coulombic model of the bound state, cannot be justified within QCD. 



Figure 6: The width T^g is shown as a continuous line as a function of T/\E±\. The cor- 
responding analytical result for T 2> I-E2I, as obtained in |12 |, is plotted as a dot-dashed 
line. 



The gluo-dissociation cross section is shown in Fig. [2j The impact of the 
colour-octet effects is dramatic for gluon momenta, q, close to the threshold: 



the full cross section (19 1 falls off exponentially while the BP cross section (151 
falls off like (|q| — ci.bpP^ 2 - At larger gluon momenta, if the thresholds in both 
cross sections are chosen to be equal to \Ei\, then for |q| = 2\Ei\ the full cross 
section is larger than the BP one by about 12% and for |q| = 5|£i| by about 
22%. In the asymptotic high-momentum limit, the full cross section overshoots 
the BP one by about 13%. The thermal decay width is shown as a function of 
the temperature in Fig. [3] Much larger temperatures than those plotted would 
likely violate the bound ma s ^ T. At temperatures lower than the energy 
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the full thermal width falls off faster than the BP one. For T = 0.3|£i|, the 
full width is smaller than the BP width by about 12%, while for T — \E\\, the 
full width is larger than the BP width by about 3% and for T = 5\Ei\ it is 
larger by about 11%. According to [13], the region \Ei\ >T> 0.3|i?i| may be 
of relevance for T(IS') produced in heavy-ion experiments at LHC. 



Note added 

While this paper was in the final writing up, a paper appeared [56] where 
the gluo-dissociation cross section in the presence of the octet potential was also 
obtained. The result is shown in Eq. (4) and for n = 1, after correcting some 
typos, agrees with our Eq. (19 1. 
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